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Abstract. In the standard model of structure formation galaxies reside in virialized dark matter haloes which extend much 
beyond the observational radius of the central system. The dark matter halo formation process is hierarchical, small systems 
collapse at high redshift and then merge together forming larger ones. In this work we study the mass assembly history of 
host haloes at different observation redshifts and the mass function of accreted satellites (haloes that merge directly on the 
main halo progenitor). We show that the satellite mass function is universal, both independent on the host halo mass and 
observation redshift. The satellite mass function also turn out to be universal once only satellites before or after the host halo 
formation redshift (time at which the main halo progenitor assembles half of its final mass) are considered. We show that the 
normalizations of these distributions are directly related to the main halo progenitor mass distributions before and after its 
formation, while their slope and the exponential high mass cut-off remain unchanged. 
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INTRODUCTION 

Understanding the structure formation process is a fundamental topic in modern cosmology. In the current ACDM con- 
cordance cosmology, the matter density of the Universe is dominated by cold dark matter (CDM), whose gravitational 
evolution gives rise to a population of virialized dark matter haloes spanning a wide range of masses [1]. Numerical 
simulations of structure formation in a CDM universe predict that these dark matter haloes contain a population of 
subhaloes, which are the remnants of haloes accreted by the host, and which are eroded by the combined effects of 
gravitational heating and tidal stripping in the potential well of the main halo. 

Understanding the evolution of the subhalo mass function, as a function of cosmology, redshift, and host halo mass, 
is of fundamental importance, with numerous applications. For one, subhaloes are believed to host satellite galaxies, 
which can thus be used as luminous tracers of the subhalo population. In particular, linking the observed abundances of 
satellite galaxies to the expected abundance of subhaloes, provides useful insights into the physics of galaxy formation 
[2, 3, 4]. Studies along these lines indicate that galaxy formation becomes extremely inefficient in low mass haloes, 
and suggest that there may well be a large population of low mass subhaloes with no optical counterpart [5, 6]. 

In principle, though, these truly 'dark' subhalos may potentially be detected via y-ray emission due to dark matter 
annihilation in their central cores [7, 8, 9, 10, 12, 1 1] or via their impact on the flux-ratio statistics of multiply-lensed 
quasars [13, 14]. Alternatively, these techniques may be used to constrain the abundance of subhaloes, which in turn 
has implications for cosmological parameters and/or the nature of dark matter. 

The evolution of the subhalo mass function is also of importance for the survival probability of disk galaxies [15, 16] 
and even has implications for direct detection experiments of dark matter. 

Finally, understanding the rate at which dark matter subhaloes lose mass has important implications for their 
dynamical friction times, and thus for the merger rates of galaxies [17, 18, 19]. 

Despite significant progress in the last years, there are still numerous issues that are insufficiently understood. What 
is the mass function of haloes accreted onto the main progenitor of a present day host halo? Does the distribution 
depend on the host halo mass and observation redshift? Answareing this question is also very important because, while 
the mass of the satellite decrease under the gravitational influence of the host due to gravitational heating and tidal 
stripping, its luminosity is thought to be preserved (author?) [20, 2 1 , 22] . In this work we address these questions using 
high resolution numerical simulations. We trace back the evolution of host haloes starting from different observation 
redshifts, zo = 0, 0.5, 1 and 2, identifing progenitor haloes directly accreted onto the main branch. 



SIMULATION AND POST PROCESSING 



Let us consider the outputs of a cosmological N-Body simulation of a ACDM universe in a periodic cube of side 
lOO/i^Mpc. This simulation has a total number of particles of 400 3 , with an individual mass of m p = 1.73 x 
10 9 /i _1 M Q . The cosmological parameters are (£2„„ a g , h, Q. b h 2 ) = (0.3, 0.9, 0.7 0.0196), more details about the GIF2 
simulation and the post processing can be found in [23] and [24]. 

At each simulation snapshot haloes have been identified with the spherical overdensity criterion. Having estimated 
the density of each particles and sorted them by density, we start from the densest one and grow a sphere of matter 
around its center, and stop when the mean density within the sphere first fall below the virial value [25]. At this point 
we assign all particles within the sphere to the newly formed halo, and remove them from the global list. We take the 
centre of the next halo at the position of the densest particle among the remaining ones and grow a second sphere. We 
continue in this manner until all particles are assigned to haloes or to the field. 

Starting from a given observation redshift zo we build the merging history tree, for all haloes 1 in the simulation more 
massive then lO 1L5 /z" 1 M , using the halo catalogues at all snapshots, separated by redshift inteval dz,, as follows. 
Considering a halo with virial mass M zo at zo we define its progenitors at the previous output z\ = Zo + dzi to be 
all haloes contributing with at least 50% of their particles to the initial system. Among them we call the main halo 
progenitor the halo providing the largest mass contribution. We repeat the procedure, starting now at z\ and considering 
the progenitor at Z2 = Z\ +dz2 of the zi-main halo progenitor, and so backward in redshift. In this merging history tree 
we term satellite all progenitors which, at any time, merge directly on the main progenitor, contributing at least with 
50% of their particles to the initial system at z = zo- 

In the hierarchical growth of the parent host haloes we identify Zf, the formation redshift of the halo, to be the 
redshift at which the main halo progenitor assembles half of its final mass [26]. Because small systems collapse first 
and merge together we expect more (less) massive haloes to have a lower (higher) formation redshift [27]. 
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FIGURE 1. Mass accreted in satellites (unevolved subhalo mass function) by the main branch of the tree at all redshifts. In each 
panel we show the results following the tree starting from different observation redshifts zo and host halo masses. The solid lines 
represent the fitting function: Equation 1. 



Haloes which exceed more then 10% the initial virial mass at any z> zo are not considered in this analysis. 



Satellite Mass Function 



In Figure 1 we show the satellite mass function at four different observation redshifts zo =0, 0.5, 1 and 2, the virial 
mass of each satellite m v is expressed in units of the host halo initial mass M Zf) . In each panel different data points and 
line types refer to various host halo masses. From the figure we notice that all data points, at each observation redshift, 
are well fitted by the following equation: 



AN 



dln(m v /M Z0 



= N Q x- a e- 62&3x \ x = 



(1) 



which is represented by the solid curve. The normalization Nq and the slope a are respectively equal to 0.21 and 0.8 
[23]. This universal distribution has been predicted by [26] using the extendend-Press & Schechter formalism and 
linear scale-free power spectrum. [28] found a different slope and normalization following haloes from zo = to high 
redshift using a Monte-Carlo merger tree realization of halos in a ACDM cosmology However the first measurement 
of the universality of this distribution in numerical simulation has been done by [23] following haloes only again from 
zo = 0. In this work we show that also if we follow haloes starting from different observation redshifts, that are zo > 0, 
the distributions still preserve their universality with identical slope and normalization. 
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FIGURE 2. Satellite mass function accreted before (left) and after (right) the host halo formation redshift. Different data points 
and line types refer to varius final host halo masses. The solid curve represent equation 1 with the correct rescaled normalizations. 

In Figure 2 we show the satellite mass function for different initial host halo masses and observation redshift, 
but considering satellite progenitors accreted before (left) and after (right) the formation redshift of the main halo 
progenitor. Equation 1 also in these cases fits quite well the data point distributions but with different normalizations 
that we term Afo.b (for z > Z/) and Afo.a (for z<Zf). 



Following the Main Branch 

In order to estimate the normalization of the distribution of the satellite mass function before and after the formation 
redshift of the host halo we need to compute the distribution of the main halo progenitor at z/- This is because, if M zo 
is the host halo mass at zo, and all the mass is assembled throught merging of progenitor haloes at zi, requiring the 
mass conservation we can write: 

Y, m m + M zi = M z > 



where we divided the sum into the total mass in satellite haloes plus the main halo progenitor. At this point we can 
also do the same considering the progenitor haloes of M Z1 at Z2 = Z2 + dz2, which gives: 



E W <>« + E m '.*2 + M Z2 = M Z0 



Following the main halo progenitor back in time until it assembles more then half of its initial mass we can write: 



E E m '^+^o = M Z : 
j=0 i 



(2) 



where jX — M Zf /M, represents the initial main progenitor mass fraction at z/ and - J^i m i,zj tne tota l mass accreted 
in satellites for z < Z/. This first term on the left hand side of Equation 2 is directly related to No M , the normalization 
of the satellite mass function in the right panels of Figure 2. 

The analytic estimate of ji can be obtained considering the distribution of the main halo progenitor mass soon after 
the formation redshift, 



with 1/2 < jj. < 1, estimated by [29] using the extended-Press & Schechter formalism. 



(3) 




0.2 0.4 0.6 



1 0.2 0.4 0.6 0.. 



FIGURE 3. Formation mass distribution measured in the GIF2 simulation for different final halo mass bins and observation 
redshifts with 1/2 < \l < 1. The various line type histograms show the result of different final host halo masses. Considering that 
the distribution does not depend on Mo, we plot all the halo more massive than lQ n 5 M (}! /h with filled circles. The corresponding 
error bars assume Poisson counts. For ji < 1/2 the mass distribution just before the formation is shown. See the main text for more 
details. 

Performing one more step along the main branch of the tree we can re-write Equation 2 and obtain the main halo 
progenitor distribution soon before its formation: 

E E m '>; +E m ^/+1 + M'A^b = M zo . 

j=0 i i 



where jj,' — M Zf+l /M zo is the initial main progenitor mass fraction assembled soon before the formation redshift. The 
jii' distribution has also been estimated by [29] and has the following equation: 

where 1/4 < ju' < 1 /2. In this case the mean value of \i' is related to the total mass accreted in satellites for z > Zf, 
which is proportional to Afrj.z,, the normalization of the satellite mass function in the left panels of Figure 2. 

In Figure 3 we show the main halo progenitor mass distribution before (1/4 < jU < 1/2) and after (1/2 < jU < 1) 
the formation redshift. Histograms with different line types refer to various initial host halo masses at corresponding 
observation redshifts zo in the GIF2 simulation, while the data points refer to all host haloes considered with mass 
larger than 10 11 5 /z _1 M Q . The two solid lines in each panel represent Equations 3 and 4 which perfectly reproduce the 
simulation results. Mass conservation implies that: 

No=No, a +N ,b, 

and for the mean main halo progenitor mass assembled soon before and after the formation redshift: 

1 =jX' ' + jx. 

The normalizations Afo.a and Afo.b have been estimated measuring the mean halo progenitor mass, for all haloes starting 
from every observation redshift, soon before and after the formation redshift: 

N , b = P-N = (0.583 ±0.09)JV 
No.b = jl'No = (0.432 ± 0.06)JVo , 

absolutely compatible within the error bar with the mass conservation equations above. The small deviation from unity 
of the sum p, + p,' is due to the smooth particle component (mass not in haloes) accreted in numerical simulations along 
the main branch. In the panels in Figure 2 the solid curves represent the satellite mass function (Equation 1) with the 
obtained rescaled normalizations that fit quite well the respective data points. 

SUMMARY REMARKS 

The standar scenario of structure formation predicts that dark matter haloes grow as consequence of repeated merges 
events. Central galaxies reside in the main halo progenitor center, while satellite galaxies are in subhaloes. These 
represent the core of surviving haloes accreted by the main halo progenitor during its cosmic evolution. In these 
pages we have measured the satellite mass function accreted by the main branch of the merger tree of various host 
haloes, starting from different observation redshift zo in a numerical N-Body simulation. The satellite mass function 
turns out to be universal: both independent of the host halo mass and redshift. The shape of the distribution remains 
unchanged also when satellite haloes accreted before and after the formation redshift are considered. In this case the 
normalizations can be obtained studying the main halo progenitor distribution soon before and after its formation. 
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